Breakdown of integrability in a quasi-one-dimensional ultracold bosonic gas 
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We demonstrate that virtual excitations ol higher radial modes in an atomic Bose gas in a tightly 
confining waveguide result in effective three-body collisions that violate integrability in this quasi- 
one-dimensional quantum system and give rise to thermalization. The estimated thermalization 
rates are consistent with recent experimental results in quasi- ID dynamics of ultracold atoms. 
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Thermalization does not occur in integrable systems 
[H , since the number of their integrals of motion equals 
exactly the number of their degrees of freedom, thus 
such a system always "remembers" its initial state in the 
course of its dynamical evolution. In an integrable sys- 
tem the finite spread of initial energy may lead only to 
relaxation towards the generalized Gibbs (or fully con- 
strained thermodynamic) ensemble Q. Strictly speak- 
ing, there is no thermalization in any closed system, but 
for non-integrable systems the eigenstate thermalization 
hypothesis 3] holds, enabling dephasing to mimic the re- 
laxation to the thermal equilibrium. 

The Lieb-Liniger model 4] of spinless bosons with con- 
tact (point- like) interaction in one-dimension (ID) is a 
prime example of such an integrable system. 

Ultracold atoms in strongly elongated traps with uo r 3> 
uj z (w r , uj z being the frequencies of the radial and longi- 
tudinal confinement, respectively) are an ideal system for 
studying ID physics as long as both the temperature T 
and chemical potential jj, are small compared to the en- 
ergy scale given by the transverse confinement: 



(1 < huj r 



krjT < flUJr 



(1) 



Strong inhibition of thermalization was observed in a 
beautiful experiment with bosons deep in the ID regime 
Q. However, recent experimental results for a weakly 
interacting Bose gas easily fulfilling the conditions of Eq. 
©ESI are in a good agreement with the thermal- 
equilibrium description of the ID atomic ensembles. 

In the present letter we investigate the breakdown of 
integrability and thermalization in ultracold ID bosons. 
The key observation is that a radially confined atomic 
gas is never perfectly ID, and radial motion can be ex- 
cited, either in reality or virtually even if Eq. {T]) holds. 
Therefore we call such systems quasi- ID. 

We start from identical bosons in a tight waveguide 
with radial frequency u r (lu z — 0), interacting via the 
pseudopotential ATrh 2 m~ 1 a s S(r — r'), where m is the 
atomic mass, and a s the s-wave scattering length: 
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We expand the atomic field operator as follows: 
W\F) = 2^ a {n,i}k ( PnA x iy) 7f — • 
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Here L is the quantization length and ip n j(x,y) is the 
normalized eigenfunction of both the radial confinement 
Hamiltonian, H^ip n j[x,y) = [n + l)luo r ip nj t. (x, y) and 
the z-projection of the orbital momentum, — i[x(d/dy) — 
y{d / dx)]ip n} i{x,y) — tf n ,t{x,y). The main quantum 
number n = 0, 1, 2, ... , and the orbital-momentum 
z-projection quantum number I is restricted by \l\ = 
mod (n, 2), mod (n, 2) + 2, . . . , n — 2, n and thus has the 
same parity as the main quantum number. The atomic 
annihilation and creation operators as n ^ j, and a\ n ^ k 
obey the standard bosonic commutation rules. 

If two colliding atoms are initially in the transverse 
ground state of the radial confinement (ID system), then 
their orbital-momentum quantum numbers after collision 
are restricted to —I and +£. 

For two-body collisions to contribute to thermaliza- 
tion, they have to lead to transverse excitations. This 
rate of population of the radially excited modes by pair- 
wise atomic collisions r2b, can be estimated for a non- 
degenerate Bose gas, using Fermi's golden rule. For 
knT < hoj r this rate is 



2^h 



n\DOt, 



- 2b ' 



2v / 2w r C e 



(5) 



where niu is the linear density and l r = \/Ti/ (muj r ) is 
the size of the transverse ground state. The dimension- 
less quantity £ = niDa 2 /l r combines two dimensionless 



parameters (niDOe s 



and a s /l r ) which can be seen 



as characterizing a ID system 

Eq. ([5]) has a transparent physical interpretation: 
is related to the 3D atomic density, which is ~ nuj/l 2 , 
times the s-wave scattering cross-section ~ a 2 , times the 
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FIG. 1: Feynman diagrams for the effective three-body pro- 
cesses in the second (a) and third (b) orders of perturbation 
theory. Solid and dashed lines correspond to atoms in the 
ground and excited states of the radial trapping Hamiltonian, 
respectively. 



exponential Boltzmann factor for the fraction of atoms 
fast enough to scatter into higher radial modes, times the 
corresponding velocity of the collision, ~ h/(ml r ). 

Calculating the numbers for the data in the 87 Rb ex- 
periments Q, a s — 5.3 nm, nm — 50 /Uin -1 , ijo r /{2ii) = 3 
kHz, T — 30 nK (£ w 0.007) one obtains a collision rate 
of T2b ~ 0.02 s _1 , at least one order of magnitude too 
small for the time scale of the experiment. 

If the kinetic energy of the collision is less than 2hui r , 
then the radial modes can be excited only virtually and 
contribute to the system dynamics only in the second 
and higher orders of perturbation theory. If after the 
collision the radial motion state is | , ^'i } , {n' 2 ,£ 2 }} = 
|{0,0}, {2p, 0}), then only one more collision is enough 
to quench the virtual excitation and return the system 
on the energy shell [Fig. [lja)]. Such a process yields an 
effective three-body collision already in the second order 
of perturbation theory. 

In contrast processes involving a virtual excitation to 
{{n'^-e}, {n' 2 ,+£}), I ^ 0, shown in Fig. con- 
tribute only in the third order, and thus will be neglected. 

The small parameter in our perturbation calcula- 
tion is ni£>a s . To avoid complications related to the 
confinement-induced resonance in ID scattering (lol | we 
assume a s <C l r Q. 

In a first step we evaluate the matrix element 



({0,0}, {2p,0}\S(x-x')6(y 



y')\{0,0}, {0,0}) = 

(2^)-i (6) 



that corresponds to two atoms in the ground state of 
the incoming channel, one atom remaining in the same 
state, but the other one being excited to a state with 
zero orbital-momentum quantum number and even main 
quantum number n — 2p, p = 0, 1, 2, ... (n and t are 
required to have the same parity). 

The result of Eq. ([6|) should not be confused with the 
matrix clement where the outgoing channel is character- 
ized by excitation of a higher radial mode of the relative 
motion of two atoms as discussed in [I~0 | , which equals to 



(27tZ^) -1 . The latter is a linear combination of the matrix 
elements corresponding to vertices on both Fig. QJa) and 
Fig.QJb), and thus can not be applied to the calculation 
of the second-order process. 

Using the matrix element Eq. ^ we can rewrite Eq. 
@ and by adiabatically eliminating the radially excited 
mode operators obtain the effective ID Hamiltonian: 
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where the summation in Eq. (jSJ) is taken over all the 
kinetic momenta obeying the conservation law k[ + k' 2 + 
k' 3 = ki + k 2 + fc 3 , a k = a {0 o} k and £ = 4 In (4/3) w 1.15. 
The relative contribution (£ — l)/£ of the virtual states 
with the excitation energy higher than 2hui r is remark- 
ably small. 

Using Eq. © we then obtain the collision rate for the 
process shown in Fig. [Ha): 



r 



da- 



^36" 



Czb U r c 



(9) 



with C; 
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5.57. Comparison to the two-body 



rates for typical experiments are given in Fig. [2l 

The result of Eq. ^ may seem counterintuitive at 
first: the collision rate is independent of temperature, 
and is proportional to £ 2 and the radial confinement u r . 

The physics behind the first observation is related to 
the fact that the collision kinetic energy is small com- 
pared to the virtual excitation energy (according to as- 
sumption Eq. |(T||). Consequently the composite matrix 
element of the second-order process should not depend 
(in leading order) on the velocities of colliding particles 
and hence on temperature (see Eq. ((6])). In addition the 
phase space volume for the scattered particles is indepen- 
dent on the incoming momenta k±, k 2 , and ^3. 

Since effective three-body elastic scattering is the dom- 
inant process the scattering rate must be proportional to 
the 3D density squared, (ni D /w 2 ) 2 . On the other hand, 
the scattering rate contains the square of the matrix el- 
ement corresponding to the diagram in Fig. 1(a), where 
each vertex is proportional to a s , therefore this rate is 
proportional to oq. The factor K/m provides the correct 
dimensionality (s _1 ). 

We can now compare the scattering rates for the two 
routes for thermalization and breakdown of integrability 
in ID systems: thermally excited two-body collisions T 2 b 
(Eq. [5]) or effective three-body collisions T^b (Eq. [9J. For 
ksT < huj r we find a simple scaling: 
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For large £ and small temperatures the scattering rate 
due to virtual excitations dominates, and can violate in- 
tegrability even when thermalization processes due to 
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FIG. 2: Ratio between the scattering rates for the two routes 
to thermalization and breakdown of integrability in ID sys- 
tems: T2b for two-body collisions leading to excited transverse 
states, and for the effctive three-body collisions. Units on 
the axes are dimensionless. £ = 0.002 (dashed curve), 0.007 
(solid curve), and 0.02 (dot-dashed curve). The points rep- 
resent the predicted ratios for various sets of experimental 
parameters from 0] (points), 0] (crosses), and (triangle). 



simple two-body collisions are frozen out. A detailed 
comparison of the two rates and their relation to exper- 
imental parameters is given in Fig. [21 The scattering 
rate due to virtual excitations of the radial modes can 
dominate over real excitations for typical parameters of 
the recent experiment 0]. 

The damping rate of an elementary excitation that 
moves in a degenerate Bose gas much faster than the 
speed of sound also has the functional form of Eq. (jSJ), 
with C^b being replaced now by Cf as t ~ 13.76. This re- 
sult suggests that thermalization is very similar also in 



12] 



the largely unexplored quantum decoherent regime 
of quasi- ID bosonic dynamics. 

To quantify the thermalization, and thereby the vio- 
lation of integrability in a quasi- ID system by the ef- 
fective interaction (|8|). we consider a non-degenerate, 
weakly-interacting (the Lieb-Liniger parameter Q 7 = 
2as/(ni_DZ^) being much less than 1) gas of bosonic atoms 
jl3| and write out the Boltzmann equation with a three- 
body collision integral taking into account the in- 
distinguishability of the particles: 



d_ 12i 2 u? r a\m 
dt h ~ V3^ 3 n 



K 



dq 



K ^K 2 - q 2 

[fko — qfk+fk- fkfk' fk") i 



dk' 



dk" X 

(11) 



k + k' + k" 



v/3 



fc± = fc + f ±^V^^ 



K = -Vfc 2 + k' 2 + k" 2 - kk' - kk" - k'k". (12) 
3 

To solve Eq. (fTTj) we use the following ansatz for the 
perturbed momentum distribution 

/*(*) = ^-exp(-k 2 /k 2 h )[l + e 4 (t)H 4 (k/k th )}, (13) 



'Ih 



where k t h = \/2mksT /h and H 4 is the Hermite poly- 
nomial of the 4th order. We choose this form, since 



FIG. 3: Dependence of the rate F[ 4 j of thermalization induced 
by effective three-body collisions on the radial trapping fre- 
quency. niD = 30 /im -1 (dashed curve), 40 /im -1 (solid 
curve), and 50 /im -1 (dot-dashed curve). 



it is the simplest nontrivial perturbation that retains 
J dk kfk — 0. Linearizing Eq. (jlip with respect to the 
perturbation amplitude e 4 {t) we obtain the exponential 
solution e 4 (t) = e 4 (0) exp(— rj^i) with the decrement 
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where the numerical constant Cu 



1.65. Tak- 



_ 64^ 

ing the perturbation proportional to a higher-order Her- 
mite polynomial H n leads only to a minor the modifi- 
cation of the numerical prefactor, leaving the functional 
dependence on the parameters of the system unchanged 
(for example, n = 5 or 6 increases the thermalization rate 
by the factor 5/4 or 13/9, respectively). Fig. [3] shows nu- 
merical values of rj^ as a function of the ID density of 
87 Rb atoms and the radial trapping frequency. 

Similarly, we calculate numerically the thermalization 
rate T 2 ^ for two-body collisions involving the real tran- 
sitions between the ground and excited radial states. 
The velocity distribution of atoms in the ground and 
excited state was perturbed in the same way, as given 
by Eq. (|13[) . the Boltzmannian distribution of over- 
all populations between the levels being kept intact. 
In the parameter range of interest we find numerically 
r 2 4 fc ] « (0.33 ± 0.03)r 2fc , i.e. 



rf| w 0.93 uj r C,e 



(15) 



The ratio of the thermalization rates for the two-body 
and three-body processes is very close to the respective 
ratio of the collision rates, shown in Fig. [2] 

It is interesting to note that we find for both processes, 
the two-body collisions to real transverse states and the 
effective three-body processes via virtual excited states, 
that thermalization in ID needs about 3 collisions which 
are able to distribute energy. This is very close to the 2.7 
collisions required for thermalization in 3D 15 1. 

For the typical parameters of an ultracold 87 Rb gas on 
an atom chip Q (uv ~ 27r x 3 kHz, n\D w 50 /im^ 1 ) 



4 



we obtain rj^ w 2 s" 1 . This thermalization rate is 
temperature-independent and much larger than the one 
calculated from the simple two-body collisions with the 



3 x 



energy sufficient to excite radial modes r?£ 
10~ 3 s^ 1 at the lowest temperatures measured (30 nK). 
The estimated rjj, is consistant with the time needed for 
evaporative cooling of a 87 Rb gas on an atom chip well 
below %u) r @,0]. 

Comparison of our thermalization rates to the experi- 
mental results of Kinoshita, Wenger and Weiss Q is more 
qualitative. In this experiment, a degenerate 87 Rb gas 
in a two-dimensional optical lattice (u> r ~ 2n x 67 kHz, 
niD ~ 10 ^m _1 ) is split by a laser Bragg pulse into two 
groups with opposite kinetic momenta, which begin to 
oscillate in a weak trapping potential in z-direction, col- 
liding each half-period of the oscillation. The experiment 
is close to the strongly interacting regime (7 ~ 1), and 
one has to take the strong suppression of three-body col- 
lisions into account. The rate of relaxation of the relative 
motion of the two groups of atoms can be written as 



(16) 



where C rm ~ C^b and rj T < 1 is the fraction of the pe- 
riod, when the two groups of atoms with opposite veloci- 
ties overlap in space. The dominant factor in Eq. (|16p is 
the suppression of three-body scattering ^(7) via atom- 
atom correlations in both initial and final states. The 
experiment [1 61 ] showed 7-fold suppression of three-body 
inelastic losses in a quasi-lD bosonic gas with 7 w 0.5 
compared to a regular BEC. However, there is a differ- 
ence between elastic scattering discussed here and that of 
Ref. [l6| where the final state corresponds to a molecule 
and the fast atom carrying away the energy and momen- 
tum released in the three-body recombination process. In 
the elastic scattering related to dissipation the suppres- 
sion enters in both the incoming and the outgoing channel 
and should be the square of the measured value. Read- 
ing the predictions for the parameters in [j| (7 w 1.4) 
from the estimates in [l6| we find a suppression factor 
^(1.4) - 4 • 10~ 4 and T rrn < 0.03 s _1 , which is con- 
sistent with the experimentally obtained lower bound to 
the relative motion damping rate > 20 s. 

The mechanism discussed here which leads to the 
breakdown of integrability in ID is to a certain extent 
similar to the virtual association of atoms to a molecular 
dimer [l7j ]. In the present discussion, virtual excitation 
of radial modes during a two-atom collision temporarily 
localize the interatomic distance on the length scale ~ l r . 
Scattering a third atom on such a transient structure of 
finite size and mass 2m leads to thermalization and vio- 



lates integrability. In Ref. 17J , collisions of a third atom 
bring "virtual" dimers, enhanced in size by a Feshbach 
resonance down to the energy shell, thus bringing about 
"quantum chemistry" in ID. 

To summarize, we identify effective three-body colli- 
sions which arise in the second order of perturbation the- 
ory as a mechanism of breakdown of integrability and, 
hence, thermalization in ID atomic gases. These pro- 
cesses are associated with virtual excitation of radial 
modes and for weakly interacting quasi-lD Bose gases 
can be dominant at ksT < huj r . Our estimations of the 
relaxation rates in are consistent with recent experimen- 
tal observations [H, @, 0] • 
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